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PREFACE 


THE aim of the present volume is to set forth a method of treating 
the geometrical optics of any given system intrinsically, that is to 
say, without introducing any artificial scaffolding and, therefore, 
without any arbitrary splitting of the entities involved in the 
problem of tracing a luminous ray through the several surfaces 
of a system. The resultant formulae, in vector language, of 
course, as the only appropriate means of intrinsic expression, will 
then be free of any unnecessary elements, and therefore, if finally 
translated into any set of scalar formulae for the practical use 
of the numerical computer, will contain no superfluous, geometrical 
or arithmetical, complications. 

Our purpose is not to treat the whole subject of geometrical 
optics, but exclusively, or almost so, that part of it which is called 
by the short name of “ ray tracing.’”’ This is notoriously the most 
laborious part of the computer’s patient work, and becomes, 
without question, a formidable task when he has to deal with 
skew rays and non-centred systems. The problem thus limited 
can be put shortly :—Given the ray incident upon any system 
of lenses, mirrors, and prisms, find the emergent ray. 

The advantages of the vectorial method of resolving it must 
not be judged by the conspicuous shortness of the resultant formulae 
alone, but also by the simplicity of their deduction, as compared 
with the usual method, and by the facility of recalling the formulae, 
or of reconstructing them if forgotten. Again, the nature of the 
proposed method is such as to make the help of drawings—which, 
especially in the case of skew rays, become, in the best standard 
treatises, exceedingly complicated—almost superfluous. This 
circumstance will be particularly welcome to readers who are not 


endowed with strong visualising powers. In connection therewith 
v 
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the often troublesome discrimination of sign will become almost 
automatic. 

As already indicated, the mathematical idiom to be employed 
to secure these advantages and to prevent the intrusion of foreign 
and artificial elements is Vector Algebra. The reader whose 
acquaintance with that natural language of space-relations is but 
slight need not, on that account, be deterred from studying the 
methods to be described. For they will involve only the rudiments 
of vector algebra: namely, the addition and the scalar multiplica- 
tion of vectors together with an occasional application of the 
vector product of two vectors. A reader who is entirely ignorant 
of the subject could, in a few hours, acquire this amount of know- 
ledge from one of the existing treatises on it. For instance, the 
first twenty pages of the author’s Veciorial Mechanics (Mac- 
millan, 1913) contain all that is required and even more. It is 
unnecessary, perhaps, to point out that the rudiments of vector 
addition and multiplication, thus readily learnt,- will be found 
useful in other connections, such as the study of electro-magnetism 
and kinematics. It is, however, permissible to urge that British 
students, above all others, should be more widely acquainted with 
the language of vectors, so powerful and yet comparatively so little 
used. For it is their own patrimony. After due acknowledgment 
has been paid to Grassman it remains true that the vector method 
originated in this country through the creative genius of the 
great Hamilton, and was shaped into a practical instrument of 
investigation by Oliver Heaviside’s pioneer work of simplification, 
seconded in-America by the independent and equally sucpcastnl 
efforts of J. Willard Gibbs. 

In order to avoid digressions in the text of thi volume it will 
be well to give here a few explanations of the symbols to be used 
and to collect the very meagre number of vector formulae which 
will be required for our purpose. Vectors will be printed in 
Clarendon, small or capital type. Thus 


Ar etc, 
will be vectors. Their sizes or absolute values will be denoted 
by the corresponding Italics, as A, ry in the above case. If r=1, 
then r is called a unit vector. The vector sum of two vectors 
will be denoted, as usual, by A+B. The scalar product of two 
vectors A, B will be written simply AB, and their vector product, 
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which will not often be needed, VAB. With these symbols the 
only formulae required will be: 

A+B=B+A; AB=BA; VAB=-—VBA, 
A(B+C)=AB+AC; VA(B+C)=VAB+VAC. 
Similarly for the product of two sums of any number of vectors. . 
If y is the angle contained between the two vectors A, B, then 


AB=AB cos y, 
and the size of the vector product VAB is 
| VAB|=AB sin y. 

The scalar autoproduct of a vector A, or its “ square,” will be 
written AA=A?— 4 
If r is a unit vector, it is evident that r?=1. The angle y between - 
two unit vectors a, b, is given by 

ab =cos y. 
The mutual perpendicularity or orthogonality of any two 
vectors A, B is equivalent to AB=o0, while their parallelism is 
expressed by VAB=o. Instead of the latter the form A ||B will 
sometimes be employed. The concurrency (or equal sense) of 
parallels and its reverse will, whenever needed, be expressly 
stated. Finally, for any A, B, C, 
VAVBC =B(CA) — C(AB). 

Scarcely more than what has here been set down will be required 
for following freely the course of the deductions. 

The main subject will, now and then, be illustrated by examples 
or exercises. These will, for the sake of better discrimination be- 
tween essentials and accessory matter, be printed in small type. 

I gladly take the opportunity of expressing my best thanks 
to Mr. Frank Twyman, Manager to Messrs. Adam Hilger, Ltd., 
fer encouraging this little work, and to my friend Prof. T. Percy 
Nunn for reading and revising the diction of the MS. My 
- thanks are also due to the Publishers for the care they have bestowed 


on the book. 
L. SILBERSTEIN. 


Lonpon, April 1918. 
RESEARCH Dept., ADAM HILGER, LTD. 
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SIMPLIFIED METHOD OF TRACING RAYS THROUGH 
ANY OPTICAL SYSTEM. 


1. The Fundamental Reflection and Refraction Laws. 


Let the unit vector n represent the local normal to the surface 
separating two homogeneous isotropic media whose refractive 
indices are » and »’. Whether n is drawn towards the first or 
towards the second medium, is immaterial, as will be seen presently. 
Let the direction of the incident ray be given by the unit vector r, 
that of the refracted ray by r’, and that of the reflected one by r’, 
the latter two being again unit vectors. The plane r, n is then 
the tmcidence plane. The angles of incidence i, of refraction 7’, 


and of reflection 2”, Fig. 1, are given, according to the very definition /o'.4 “0: 
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of the scalar and the vector product, by as 
’ hj Cs / 
——- cosi=m, cosi’=r’n, cosi”=r’'n, a. os 


or by the sizes or absolute values of the corresponding vector | 


products, viz. sini=|Vrn|,. # P17 (46) bang. 


and similarly for the dashed angles. Thus the law of reflection, 
which asserts the equality of z, 2” and the coplanarity of r, n, r”, 
can be written, regard being paid to the sense of r”, 

Vr"n = Vin, oe RE Reg 


* The second equation, besides being obvious, is also a consequence of the first. 
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Multiply both sides of the first equation vectorially by nu. Then, 
by the last formula of the Preface, and remembering that n?=1, 
the vector form of the law of reflection will become 


r” =r — 2n(nr), (I) 


giving the reflected ray r” in terms of the incident ray r and the 


incidence normal n. 
Again, the refraction law, which asserts that p’sin7z’=psint aa 


Evia that r, n, r’ are coplanar, is tantamount to eee 
Lge Vns’ = Vns, foe ee 2 oli he. Ja 


Pee 
‘where s and s’ stand for wr and p’'r’, respectively, And might, ae 4 


therefore, be called directed light paths. The name is of small 
importance. Keep in mind, however, that the r’s being unit 
vectors, 8?=y? and s’=”. Treat the last equation exactly as — 
in the case of reflection, z.e. multiply it vectorially byn. Then the - 
result will be s’=s +n(s‘n—sn), 
b F a a i es 
or, putting g=sn-—sn=p’' cost’ —p cost \Fee +> fF / 
as a convenient abbreviation, to be used throughout the book, 
s’=s+gn, (II) 


This is the required vector form of the law of refraction,* giving 
the refracted ray s’ in terms of the incident s and the normal n. 
It is true that the scalar, or ordinary numerical, coefficient g itself 
contains the unknown vector s’, but only in the combination 
sn=p’'cosz’. Now n,r being given, 7 is given; whence 7’, and 
therefore also g, can be rapidly calculated by means of the relation 
pe’ sini’=psinz. It is, therefore, not advisable to eliminate 
s’ formally from the right-hand member of (II), although this 
can easily be done. (Vide infra, Section on Multiple Reflection.) 
It will be found more convenient to leave the form (II) of the 
refraction law with the numerical coefficient g as it stands. 


* Formula (II) has also been arrived at, in a roundabout way, by Sommerfeld 
and Runge, Anz. d. Physth, 35) I9II, p. 277, the object of the authors being to 
give a general differential ‘basis to geometrical optics. In our case it appears 
simply as the immediate vector translation of the law of refraction considered as a 
given law. Cartesian (three formulae) equivalents of (II) are, of course, to be 
found also in some of the old works. Formula (I) and its consequences were 
given, in dyadic form, in my paper ‘On Multiple Reflexion,” Phil. Mag, for ~- 
November 1916, 
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A glance at (I) will suffice to see that the law of reflection follows 

from that of refraction, (II), by putting in g 

—p instead of p’, (III) 
a well-known rule. In fact, it is enough for this purpose. to multiply 
(1) on both sides by the refractive index u of the medium common 
torandr”’. Remembering the simple rule (III), we can henceforth 
use (II) alone as the more general formula. Whenever a separating 
surface is contemplated as a reflector write g= —2sn= — 2m cos 7. 
Unléss the contrary is expressly stated we shall conceive every 
such surface as a refracting one. 

Notice that if in (II), and therefore also in (I), n is replaced by 
—n, the formulae remain intact. 

Let now the refracted ray impinge upon another refracting 
surface. Then the new refracted ray will be obtained from the 
tay s’ as this was from s, and so on, for any number of surfaces. 
Call n,, ny, etc., in general n,, the unit incidence normals of the 
consecutive surfaces (at the places of incidence of the ray in 
question) and pe, u,’ the indices of the media separated by the 
«th surface (Fig. 2), so that py’ =Pe41- 

/ 
By Bb; =f. /2=h3 etc. 


r 
Fic. 2. 


Further write, for each surface, 
Be=Pidey, 8 = fe Tes 
as above for a single surface, and put g.=8,.—82,. Then, 
applying the law (II) to the successive surfaces and remembering 
that s,’=8,+41 (newly refracted = actually incident ray), and writing 
s,=s for the original incident ray, the finally emergent ray will be 
given by s’=s+ 20, | (IV) 
the sum to be extended over all surfaces. 

In the case when all the surfaces are reflectors the last equation 
assumes a remarkably simple form, as will be seen irra later section, 
and offers as a final formula conspicuous advantages. In the 
case, however, of refracting or mixed systems (refracto-reflectors) 
we shall not actually require the resultant (IV), but shall employ 
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the formula (II) successively at each surface, after having developed 
an equally convenient formula for the ‘‘transfer’’ of the data 
concerning the last refracted ray to the next surface. This latter, 
the transfer formula (which will replace and supersede * the formid- 
able array of what the German opticians call ‘“‘ die Uebergangs- 
formeln’), is an unavoidable supplement of (II), the refraction 
formula. The former will be soon developed in one of the fol- 
lowing sections, again without the “‘help”’ (which is rather a 
nuisance or obstacle) of any scaffolding or framework of reference. 


2. The Nature of the Data defining an incident Ray. 


It will be recalled from the rudiments of vector algebra that all 
of our vectors, such as r or 8, are free or non-localized vectors. 
That is to say, any line segment of length s and concurrently 
parallel to a vector s drawn concretely somewhere, is, by the very 
definition of the equality of two vectors, for all purposes of our 
vector algebra, equal tos. In other words, to gives means: to give 
#and the direction of the ray only. Now, if the refracting surface 
is a plane, this is all that is needed to find the refracted ray s° 
(u’, of course, being given). For the whole plane t has but one 
normal n. Thus, unless we desire to know, for instance, the glass 
thickness traversed by a ray in a prism, the localization of the 
incident ray s is immaterial and not needed. Notso, however, when 
the refracting interface of the two media is a curved, say a spherical, 
surface. For then n is a function of position on the surface. Thus, 
s giving the direction of the incident ray, we require further data 
for its localization. Now such a datum is obviously the place 
on the surface where it is struck by the rays. This amounts to two 
scalar data defining position on the surface, let us say uw, v, measured 
along any network of curvilinear coordinate axes spread out over 
the surface, as for instance, its lines of curvature. We need by 
no means specify these u, v. Clearly, ‘‘ to give the incident ray’ 
means to give s{ and the place (u, v) of its incidence upon the 

* In the sense that it will be seen to be applicable to cases to which the usual 
formulae are not. 

+ Euclidean space being, of course, tacitly assumed throughout. 


{Any ¢wo scalar data, the setter of the problem prefers. We do not count » 
involved in s=yr. In general we shall not count data describing the optical 
system, vzz. refractive indices, curvatures, and so on. 
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surface in question. But this being given, the normal n is also 
‘given, which—being a unit vector—amounts precisely to two 
scalar data. 

Thus, whatever the nature of the numerical data for the incident 
tay, be they angles or projections upon any lines or planes, or 
distances from axes (artificial or intrinsic), we are fully entitled 
to say that “incident ray given ” means, for a surface o,, 


S, and 0, given. 


Such being the case, formula (II) finds the refracted ray s,’. 
Moreover, it finds it together with its localization; for the latter 
ray starts from the point of the surface o, struck by the incident ray. 
But in order to apply (II) to the next surface o,41 we must transform 
the last acquired data s/=8.i1, 0 

into Seftyo Nes1; 

that is to say, we must learn how to replace n, by nx41. It is this 
which will be given by the transfer formula. Manifestly, having 
once obtained 8,41, 0x41, we stand before the surface o,+1 precisely 
in the same position and with the same equipment as before the 
_ preceding surface. Thus the process of tracing the ray can be 
- continued indefinitely. 

There is no essential difficulty in developing, in correspondingly 
general terms, the transfer formula for a pair of surfaces of any 
form. But for all purposes actually met with in technical optics, 
it will be enough to develop the formula for spherical surfaces (of 
any finite radii), and later on (in a special section on prismatic 
systems) also for the sub-case of two planes or a plane and a sphere. 
In the more important case of spheres proper the transfer formula 
is indispensable for finding even only the direction of the refracted 
ray, 8’.+1, While in the case of planes (prisms) the ray can be traced 
through the whole system, as far as its direction is concerned, 
without the “transfer formula”’; the latter being required only 
when one is interested in the thickness of the several media 
traversed. 


3, The Transfer Formula for Spherical Surfaces. 


‘It was said in what precedes that, the point (wu, v) on a surface 
being given, the incidence normal n is thereby given. In general, 
however, the converse manifestly does not hold; two or more or, 
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as in the case of a plane surface, even ? points having parallel 
normals. But if the refracting surface is a sphere, of centre O, then 
we have a one-to-one correspondence of points (u, v) and normals n, 
there being no two points to which the same n belongs. This 
offers the advantageous possibility of considering the incident 
ray to be given directly by s and n. It will be convenient in the 
case of any spherical surface to draw all the normals n as unit 
vectors from the sphere’s centre as their common origin. The 
two vectors s, n being given, we know already, by (II), how to 
trace the ray through the corresponding surface, i.e. how to find 
the refracted ray s’. 

What we still require is to make this ray ready for tracing it 
through the next surface of the system. This is an easy geometrical 
problem, although the optical authorities have made it a compli- 
cated one by treating it with too many “ auxiliaries.’”’ Treating 
it (according to our principle) intrinsically, it is, and will remain, 
a simple problem. 

Let oj, v,, as representatives of any pair ox, 0.41, be two con- 
secutive refracting spherical surfaces of the given optical system. 
Let O,, O, be their centres, and R,, R, their radii. Draw all the 
unit normals n, from Oy, and all n, from O,, asarule. The optical 
system, as always, being given, so are O,, O,, and all the following 
centres, Og, etc., collinear or not (centred or not centred system). 
Thus let D be the scalar distance of the two centres contemplated 
and a an unit vector drawn from O, towards O,, as in Fig. 3. (The 


Fic. 3. 


reader should note that the surfaces need not turn their convex 
sides towards one another, nor need QO, precede 0,; the formulae ° 
to be given presently will hold always.) Let n, be the particular 
normal pointing to P,, the place struck by the incident ray s, 
which penetrates into the second medium as 8,’ =sg, and strikes the 
next surface in P2, a point to which the particular normal-n, 
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corresponds. Having found s,’=s, by (II), we wish to find P,, 
t.é. the normal n,. 

In order to obtain it walk round the quadrilateral 0,P,P,0, 
back to the starting point 0,. This being a closed polygon, no 
matter whether plane or not, or whether you recross your path * 
or not, the vector sum of the four vectors traversed is nil, 7.e., 
calling L the (scalar) length P,P,, or medium “‘ thickness ” traversed 


by the ray, and remembering that s,’=s,=sro, and therefore, 
— 


Ron,=Ryn,-Da+Lr, * (1) ¢ 


This, in fact, is already the required transfer formula, expressing 
the unknown n, in terms of things given—all but the scalar co- 
efficient L, a length. In order to find this (which, by the way, is 
not only needed for “‘ tracing ” the ray but will serve the computer 
in his further task of calculating “‘ astigmatism and curvature ”’), 
eliminate from (1) the principal unknown n, by squaring both sides 
of this vector equation. In doing so remember that n,?=1, as 
are also the squares of the remaining three Clarendons. Then the 
result will be 


L? —2(Pr,)L =R,? — P?, (2) 
a quadratic equation for L, with 


which is a given vector. This finds L and completes the required 
transfer formula (1) which can be written, more shortly, 
R yn,=Lr,—-P. (1’) 
Is not this simple? In order to estimate the simplification of 
the deduction, due to the vector method, read the same subject, 
for instance, in the well-known collective work of the scientific 
Jena staff,t where it covers ten pages of print (in which, to be 
just, not one but three alternative ‘‘ Uebergangs’’-methods are 
described, none inferior to the others either in rigour or in clumsi- 
ness). In harmony with this, their “ refraction-formulae”’ are 
equally complicated, the whole apparatus consisting of something 
like fifteen or sixteen formulae. The complication in the scalar 


* Which would happen with the choice of another exemplifying figure. 


+ Die Bilderzeugung in optischen Instrumenten, edited by M. von Rohr, Berlin, 
J. Springer, 1904, pp. 52-62. It is, nevertheless, a valuable work. 
B 


BT) Ar. : p23 aa / 
P,P,=Lr,, th FP eA om b% tet. 
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treatment is of course due to the circumstance that skew rays 
are aimed at, not cutting 0,0,. But as the reader will have 
noticed, we have nowhere assumed that the incident ray has any 
such privileged orientation. To repeat, the quadrilateral O,P,P.0, 
need not be plane; the resulting formulae (1) or (1'), and (2) 
with (3) are perfectly general. 

The solution of the quadratic equation (2) will, in numerical 
cases, cost but little trouble, and this will be repaid by obtaining 
L not only as an auxiliary for the transfer formula, but also as 
one of the “elements” needed for estimating astigmatic and 
associated defects of the optical system in question. The equation 
(2) has, of course, in general two roots, 

L =(Pr,) + /(Pr.)?+R.?—P? (4) 
corresponding to the two points, P, and P,’ (of Fig. 3), in which 
the ray cuts the second sphere, and in practice there will be no 
difficulty in discriminating between the root or point (P, in the 
case of Fig. 3) to be retained and that to be rejected. And if 
both roots coincide or are complex (conjugate), there is an end 


to the ray-tracing; for in these cases the ray does not penetrate 
into the next medium. 


-It would be entirely superfluous to enter at this stage into the 
details depending on the configuration of O,, O,, and the following 
centres, together with the sense of the curvatures (concave, convex ; © 
convex, convex, andsoon). The reader will do well to draw several 
figures for the possible cases, especially for those in which the sides 
of the quadrilateral cross one another,—and to inspect them care- 
fully. Among other things he will then convince himself that, with 
the vector a drawn always from O, to O, (whether O, precedes O, or 
follows it), formula (1), supplemented by (2), (3), remains always 
valid. 


Having obtained n, by means of the above formulae, we can 
proceed to trace the ray through the surface having its centre 
at O,; similarly we shall pass from that to the surface centred 
at Os, and so forth. No complications will be met with if it should 
happen that the system is not centred, i.e. if 01, 02, Os, O4, etc., 
are not collinear, or even not coplanar, For the first two surfaces, 
@ =A, will be the ‘‘ axial”’ unit vector; for the second and third. 
another unit vector ag, drawn from O, towards O3, will replace 
it, and so on. These vectors, being not artificial but intrinsic, 
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inherent in the optical system, cannot give rise to any superfluous 
complications of the calculator’s work. If the centres happen 
to be arrayed in a straight line, all the better: in that case a single 
unit vector a will perform the function of the axis throughout 
the tracing process. The calculations will become correspondingly 
simpler, again for an intrinsic reason, to wit, because the optical 
system as such is a particularly simple sub-case of the general 
non-centred system. 


4. Summary of the Tracing Process, for Systems of 
Spherical Surfaces. 


It will be well to summarize the whole procedure of tracing a ray 
through any such system. 
The incident ray being given in direction and position, the data 
of the problem are, in whatever form, 
8, and ny. 
From these, 4, cos 7;=8,n, and 2,’ by p,’ sin 7,’=y, sin7, follow 
immediately, whence the scalar coefficient 
£1=/1' COS 2,’ — py COS 1,.* 
This being found we have, by the refraction formula (II), 
8; =8, +20, 
which is also s,, the ray incident upon the second surface of the 
system. 
The vector drawn from O, to the next centre O, being Dy, = D944, 


and 
P=D,,— Rin, 


the medium thickness L is calculated as the appropriate root 
L=(Pr,) + /(Prq)? + Ry? — P? 
of the quadratic equation (2), with r, =89/p2. 
Finally, by the transfer formula (1) or (1’), 
I 


Thus the two vectors are found, 
8 and Nay, 
ready to trace the ray, thus localized, through the second and 


* Which is to be replaced by — 2, cos7z,, when the surface is contemplated as a 
reflector. 
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to carry it over to the third surface, and so forth, up to the last 
surface of the system, giving the finally emergent ray through 
two vectors, say, : 


‘ 


s’ and n’, 
the first determining its direction, and the second localizing it 
completely with respect to the system. 

The results, at any intermediate surface o,, as well as at the 
last surface of the system, are all of the same form as the original 
data (s,, 2,). s 

Each of these pairs, whether of data or results, amounts to the 
knowledge of four scalar numbers, two involved in s (not counting 
#4) and two in the unit vector n. 

This exhausts the matter in its general aspect, the procedure 
being valid for any orientation and position of the incident ray 
and for any system of spherical surfaces, centred or not. 

With regard to its particular or practical aspect, the successful 
application of this compact set of vector formulae to numerical 
cases will depend only upon one negative condition, which is: not 
to spoil their intrinsic simplicity by inappropriate splitting or 
decomposition of the vectors. A large portion of the remainder 
of the present work will be dedicated to an attempt to show the 
reader how to avoid this danger. 

As a matter of fact, readers skilled in the handling of vectors 
will, after what has been said, hardly require any further explana- 
tions. This little book, however, being chiefly intended for those 
who are not experts in the use of vectors, it has seemed indispensable 
to add some hints at least which are likely to help the majority 
of readers in applying the general formulae to more or iess concrete 
cases. These will be given in some of the following sections. 


5. Hints and Examples. 


Having hitherto, in developing the fundamental and general 
part of the subject, followed a systematic course, we can henceforth _ 
without danger of confusion depart therefrom and choose our 
subjects freely. Our aim being to explain the general procedure 
by illustrations and hints, a pedantic observance of order would 
only hamper us unnecessarily. 

Let us start with the simple case of a lens consisting of a Rt 
piece of glass of index p limited by two spherical surfaces c,, o,, 


HINTS AND EXAMPLES. ax 


the surrounding medium being air or empty space. To fix our 
ideas let us take a 

Biconcave Lens. Let, as before, 0,, R, be the centre and the 
tadius of the first surface, upon which the incident ray s=r 
impinges, and O,, R, the centre and the radius of o,, the back 
surface; finally let n,,n, be the unit normals, drawn from the 
centres towards the points at which the ray enters and emerges, 
as in Fig. 4, which however is scarcely needed. By assumption, 


Fic. 4. 


the two unit vectors r,m are given, amounting each to two 
numerical data. The first praetical question is: How are they 
given? Generally speaking, this will depend upon the decision 
of the person who sets the problem. Yet, since we are here setting 
our own problems, let us put the data in a reasonable form, seeking 
only to avoid artificial frameworks of reference. Now, the lens, 
consisting but of two spherical surfaces, offers us but one intrinsic 
directed element, v7z. its axis, which we will represent by the unit 
vector a drawn from O, towards O,,—these points themselves 
having already been exploited by taking them as the origins of 
the normals. Besides a there is nothing intrinsic to build with 
a framework.* Now, if the ray r happens to lie in a meridian 
plane (passing through the axis), then nothing more than a is 
needed for the ray’s relevant definition; say the angles r and n, 
make with the axis. But let r be a skew ray, not contained in any 
meridian plane. Then the lens itself offers us manifestly nothing 
to complete the determination of the direction r. But the ray 


* A plane normal to a may be a welcome auxiliary, but gives, manifestly, nothing 
independent of a itself. Nor shall we look for the rim of the lens to help us. 

+ The simpler case of a meridian or non-skew ray may be taken up as a sub-case, 
and as an easy but useful exercise. 
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itself, as localized, does offer us something, and exactly as much 
asis needed. In fact, it strikes the first surface in a given point Pj, — 
given in one way or another. Through that point draw your 
‘first meridian’ plane, marking it with any unit vector b, say 
perpendicular to a, and count, from that plane or semiplane, the 
geographical longitude ¢, the companion of the polar distance 
or colatitude 6,=<+a,n,. Thus 6, will be a polar system of 
coordinates with centre O,. Remember that we have drawn our 
first meridian plane a, O,, b, only through the end point P, of r, 
but that the ray is in general not contained in that plane. 

Now, this is all we want. The vector a being furnished by the 
lens, and b by the incident ray, we have everything which is- 
required to write down numerically (scalarly) the data. In fact, 
‘‘y given ’’ means now ra, rb given, say : 


ra=cosa, rb=cosf. (5) 


For symmetry, although it is superfluous, introduce a third unit 
vector c= Vab (normal to a, b, towards the reader), and write 


te=cosy, so that (ra)?+ (rb)? +(re)? =1. - (6) 
Keep in mind, however, that the first two are by themselves suffi- 


cient. This settles the question of r. Next, ‘“‘n, given” now 
means 0,, $, given, 7.e., to be general, 


n,=acos 6, +(bcos ¢,+¢sin $j) sin ,, (7) 
and since, by the choice of the first meridian, ¢, =0, 
n,=a cos 0,+bsin 0;. (7a) 


In fine, the complete numerical data, immediately attackable 
by the computer, are 


61, $,=0, ra, rb, 


the latter two as in (5), or, if he prefers, with re as in (6) instead 
of rb. 


Such being the scalar data, act according to the compact pro- 
cedure of Section 4. Thus * 


COS 1) =F4N, = Cos 6, (r,a) +sin 4, (r,b) 
=cos #;.cosa+sin 6, cos, by (7a), (5); 


x 


SREY tral > aes : : 
whence sin 2,’ — sint; and g,=pcos1,'—cost,, 


* Remembering that in the present case 8)=ry=Y, =I, py! =Hg=My fg’ =I. 
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ready for the first refracted ray 

pr)’ =1, +g,0y = pT. (8) 
If the reader wishes to split this at the present stage he can do so 
at once, multiplying (8) scalarly by a, bore, But, being interested 
only in the finally emergent ray, let us proceed according to the 
general prescription of Section 4. To ‘transfer’? the ray fp, 
found by (8), to the surface o,, write D for the distance 0,0,, and 


P=Da-—R,n,, as before; then the scalar product required for L 


will be 
Pr, = D(r,a) — R,(r.n,) 


D R : 
= 7, (C08 a+ g, C08 6) —— (cos 4 +81), by (8), (5), (74), 
and, the square of P, 
P? =D? + R,2-2DR, cos 6,, by (74). 

This finds the glass thickness traversed, which in the present 

case is, manifestly, the smaller of the two roots of (2), 
L=(Pr,) —V(Pr,)? — P?+R,?, 

which the reader will properly contract by using the last two 
expressions. The ‘‘transfer’’ is thus completed, giving n,, and 
therefore also the point P., where the ray strikes the back surface, 
I 

To know the numerical coordinates of this point may be of 
immediate interest. Now, calling 6, $, its colatitude and its 
longitude, with O, as centre, but still using the previous first 
meridian plane, in short the same unit vector b, we have just, 
as in (7), 


y= 


n,= —acos 6,+ (bcos $2+¢ sin #4) sin 9,, 


1.e. for the required two angular coordinates, 


cos 0,= —n,a, cos,sin#,=n,b, sin dsin 6.=n,¢, 
or, ultimately, cn, 
cos 0,= —an,; tan ¢,=—%. (9) 
bn, 


Notice that ¢, being zero, ¢, need not necessarily be so; in 
other words, n, is in general not coplanar with a, ny. It may be 
well to recall here that, for any three vectors, 


AVBC =BVCA =CVAB 
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is the volume of the parallelepipedon constructed on A, B,-C as 
edges, and that, therefore, the sufficient and necessary condition 
of coplanarity of n,, a, n, is the vanishing of the (scalar) expression 
R,aVn,n, = LaVn,r,—aVu,P, 
or, since P is itself coplanar with a, n,, of the expression 
paVn,r,=aVn,r,, by (8)- 

We thus see that ¢,#0, i.e. that P, is outside the meridian plane, 
unless the incident ray r, happens to be in a meridian plane. For 
any skew incident ray we have ¢,#0. 

The coordinates $9, 9, of the leaving point P, being thus found 
by (9), the emergent ray is localized, and it remains only to find 
its directionr,’. This is, by the refraction formula, 


\ To =MTgt+ZoMs, f2=COS%y’ —# COS 15, (10) 
the latter to be found at once from the known 


I 
TeNo re A —Pr,), 


as was g, a moment ago. This completes the tracing of the ray. 
For, if a’, B’ be the angles which the finally emergent ray makes 
with a, b, we have only to multiply (10) scalarly by a, b, to obtain 
cosa’=r,'a, cos A’ =r,'b. 
Thus the original incidence data having been put in terms of 


A, $1=0, 4, B, 
the final answer is obtained in similar language, viz. 


Oe, fo, a’, B’. 

This essentially exhausts the problem of a simple lens. If. 
instead of a biconcave we have a biconvex, or a ‘‘ vex-cave ’’ lens, 
etc., the only modifications needed will concern the sign of a 
in the expansion of ng, and the sign of the square root appearing 
in the ‘‘ glass thickness’ L traversed. The exhaustive discrimina- 
' tion of all possible cases, which will among other things depend also 
upon the relation of R,, Ry, D, may be left as an easy, but useful, 
exercise for the reader. The central glass thickness is, in the 
above case, D—(R,+R,), and similarly, mutatis mutandis, for the 
remaining cases. 

Any further information concerning the emergent ray may, 
whenever required, be easily deducted from the obtained 05, $.; a’, B’, 
or vectorially, from ng, r,’._ Thus, for instance, to decide whether 
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the emergent ray meets at all the axis of the lens (at a finite or 
infinite distance), test the triple product 

p=—aVry Dg, (11) 
which is a scalar, or ordinary number. If p=o, the emergent 
tay or its prolongation meets the axis somewhere or is parallel 
toit; and if p+0, thenr,’ isaskew ray. By (10) and (8) we have 


p=aVn,(r, +g,n,), and so on. 


Let the reader complete the development of p, and discuss the 
circumstances of its vanishing or not. A useful exercise. 


As a further exercise, assume (abstractly or in a numerical example) 
that p + 0, so that the emergent ray is a skew ray, and find, vectorially, 
by using r,’, ng, the shortest distance of the emergent ray from the 
lens axis; also the point of the axis which lies at that distance 
from the ray. 

As another exercise treat the trivial case of =o and find, again 
vectorially, the distance, from the back surface, of the point of 
intersection of the ray r,’ and the axis. 

Construct a handy formula for the angle between the original 
incident ray and the finally emergent one, remembering that its 
cosine is r,r,’. ; 

The reader will also find it useful to invent for himself similar 
problems, resolve them vectorially and evaluate the appropriately 
scalarized expressions numerically. In this way only will he acquire 
the needed skill in handling not only vectorial optics but vector 
investigations in general, which, after all, will be found also to be a 
comparatively pleasant occupation. 

Not confining himself to lenses, let him also treat a single or two 
spherical mirrors, for which g, assumes the simpler form —2cosi,, 
as explained before. 

In actual numerical ray tracings the principal and the accessorial 
vector formulae will best be written out (and glanced at whenever 
necessary) in the margin of the computer’s sheet, or on an independent 
little sheet of paper conveniently placed before his eyes. Now and 
then it will be found convenient to note a vector result by writing 
out the figures found for its ‘‘ components,” as for example 

P=15a—24-75N1, 
having previously decided upon the choice of the length unit. It 
will be remembered, of course, that every wnit vector is dimensionless, 
being a vector divided by its own (scalar) size or representative 
length. Finally, details concerning the manner of writing and 
arranging the numbers or their logarithms will best be left to the 
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individuality of each reader. It is a familiar experience that, within 
very wide limits, everybody finds his best way of arranging the 
numbers by himself, after he has repeatedly executed a given type. 
of computation. It is by no means desirable to drill away human 
individuality. 


6. Prismatic Systems. 


Leaving for the present the spherical separation surfaces, let us 
turn to plane ones. A series of different optical media separated 
from one another by plane surfaces o, constitutes a prism—or 
prismatic system. This is of almost equal interest and importance 
with lens systems. As a matter of fact, any prism system is but 
the simplest sub-case of what has been treated in the first two 
sections. Yet, in order not to omit any opportunity of familiarizing 
the reader with the vector treatment of optical problems, let us enter 
somewhat into the details of these comparatively simple systems. 

A single Prism. Let us begin with the case of but two plane 
surfaces o,, 7, limiting (to an extent) a medium, say glass, of 
refractive index ». Outside 4=I1. In the first place, let us treat 
the rays apart from their localization, i.e. only as far as their 
directions are concerned. Each refracting face is then fully 
characterized by a single unit normal, n, for the first, n, for the 
second plane face. Remembering that inversion of the n’s leaves 
all the fundamental formulae intact, draw both normals towards 
the air, as in Fig. 5, so that the prism angle or refracting angle e 
is given by 

ny, =cos(m —«) = —cose. (12) 

The question of an ‘‘ auxiliary” framework now entirely dis- 
appears. Unlike the simple lens (or any centred lens system) the 
prism is liberal enough to furnish us, by its own intrinsic elements, 
a complete natural reference system. For it has, together with n,, ng, 
a well-determined edge. Take a unit vector e along the edge, 
normally to the plane of Fig. 5, and towards the page (or down- 
wards if your reading table is horizontal), so that 


Vn,n, =e sin (7 —e) =e sine. (13) 


The angle ¢ being given, of course, this is a simple relation between 
the three natural reference vectors. Returning to Fig. 5, drop 
the outline of the prism itself, shift the normals n,, n, to a common 
origin, and draw instead the unit edge e, as in Fig. 5a. The prism 
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will then be represented by the simple three-vector skeleton, with 
p. subintended. 


Now let any ray r impinge upon (the “‘ first ” face of) the prism 
thus represented. This given ray may also be drawn from the 
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same origin, though it is scarcely necessary. It finally emerges 
from the prism as r,’=r’, say, viz. by (IV), Section 1, 


r’=r+gy0j + goNo. (14) 
Here r is, of course, in general any skew ray. The explicit 
expressions for the scalar coefficients g,, g. can be found at once. 
But even without the knowledge of their values the vector formula 
(14) shows us at a glance an interesting and perfectly general 
property of the prism. Multiply it scalarly by e; then the result 
will be 
r’e=re =cos @, say. (15) 
Thus the final emergent ray makes with the edge the same angle 0 
as the original incident ray, whatever its orientation. In other 
words, the only effect of the prism is to turn the ray rigidly round 
the prism edge as axis. The amount of this rotation depends upon 
the properties of the prism, » and ¢«, and upon the orientation of 
the incident ray. 
The ray r being given, so is the incidence angle 7=arccos (rn)), 
whence also the coefficient g,, appearing in the resultant formula 


(14), 
g, =p? — sin? 7 — cos 7. (16) 


Similarly the second coefficient g, could be written down, as in the 
case of the lens (and somewhat easier). But it is more elegant to 
represent fully the final emergent ray r’ without the help of g,.* 


* Especially as the orderly way of calculating g2 would presuppose the know- 
ledge of the final exit angle 2’=7,'. 
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To get rid of this troublesome coefficient multiply (14) vectorially 
by n,; then, by (13), and remembering that Vn,n,=0 identically, 


1 &1 
Vr'n, = =Vm,+5"° : (17) 
This expresses a vectorial property of the finally emergent ray 
in terms of things which are all known. Together with the capital 
property (15) the last equation completely determines the emergent 
ray in its orientation. y 


Exercises. Discuss equation (17), deducing from it various 
vectorial consequences and translating each of these into geometrical 
language. 

Find the expression for the amount of rotation round prism edge 
as axis, which converts the incident r into the emergent r’. 

Find, in the general case of any skew r, the limit condition, beyond 
which the ray does not emerge but is internally reflected (at o,). 

Treat vectorially the trivial, but instructive, caseofrie. (Assume 
at the outset r of the form an,+bn,; a, 6 scalars.) Translate your 
results into scalar language, deriving thus the usual formulae. 

In practice (spectroscopy, etc.) the aim is, of course, to make 
re=o, or rile. Yet the manufacturer, and the student in his 
physical laboratory, can but approach this ideal within, say, +10” 
or even broader limits. It is, therefore, important to estimate 
the effect of this error upon the emergent r’. Develop, therefore, 
the above formulae up to the very end by assuming a value of 


re=cos 0 
slightly differing from zero, 1.e. equal to a small fraction. In other 


words, put 9== +3, where 6 is a small angle, say 10” or 1’ or 5’, 


Employ, whenever necessary, the well-known developments in series : 


xs 
sin #=%— st us 
cos #= 1-5 +5- 


Find the effect upon r’ of a small change or “ variation ”’ in the 
direction of the incident ray r. Replace r by r+6r, remembering 
that r?=1, and therefore, r r=o, or dr |_r, and so on. 

Similarly study the effect of a slight variation Se in the refracting 
angle of the prism, as expressed by (12). Scalarize the resulting 
formulae and interpret them in physical (geometrical) language. | 


Medium Thickness traversed. In what precedes, the rays 
incident upon, transient through and emerging from the ptism 
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were considered with regard to their direction only, independently 
of their localization. In fact, it was possible for that purpose to 
replace the prism by a mere triad of unit vectors, n,, mg, e. The 
“transfer formula’? was not needed owing to the circumstance 
that the same ny, served as incidence normal for the whole second 
plane. This is the reason why the treatment of a prism is con- 
spicuously simpler than that of a lens. In short, the tracing of ~ 
rays through a prism as far as their direction goes requires no 
“transfer” process. Some device of that kind, however, becomes 
ory if we desire to find the length of the ray, whose direction 
is r;’=r,, within the prism, or, as we will call it for brevity, the 
medium thickness traversed. Let this be L. 


Fic. 6. 


Obviously, in order to find it, we must localize the incident ray r 
with respect to the prism, by giving, say, the point P, where it 
strikes the first face 7,. Since P, cannot be defined by n, (as in 
the case of a spherical interface), the normal being common to all 
points of the plane, it is indispensable to determine position on oy 
by some coordinates. As such we may conveniently take 4, 
measured along e (the edge), from an arbitrary origin O on the 
edge, and x, perpendicularly to the edge. Similarly y,, x, for the 
second face o, of the prism. We can collect these pairs of co- 
ordinates by introducing the vectors 


Py =%k, +710, Po=%K2 +2, (19) 
where k,, k, are two unit vectors along the axes of 4, %2, as in 


Fig. 6. 
To obtain L, walk round the triangle OP,P,0 ; then the result 


will be 
P2=P, + Lry. (20) 
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This, in the present case, is the “ transfer formula,” somewhat 
simpler than in the case of spherical surfaces. Multiply (20) scalarly 
by n,, which is normal to p,; then 


Lyn, = — PDs, (21) 


which finds L, and completes at the same time the prismatic 
transfer formula. Also, of course, L? =(p, —D,)?. 


. Exercises. In (21) substitute r,=r,’, using the fundamental 
refraction formula, thus completing the little investigation. 

Study L in its dependence upon the prism angle «, and upon the 
direction and the position of the original incident ray. 

Deduce the formulae for the simplest sub-case re=o. 

Find the deviation of the ray, noting that its cosineisrr’. Establish 
the conditions of minimum deviation. 


System of any Number of Prisms. The tracing of any ray, 
in direction and position, by the successive application of the 


refraction formula 
Sc+1 =S« =S8_ + Lele 


and of the transfer formula typified by (20) with (21), offers no 
difficulties, nor essentially novel points, and need not, therefore, 
detain us here. We can imagine that by the application of the said 
process all the scalar coefficients g, are already determined by the 
given incident ray r and by the properties of the given system 
of prisms. 

Thus, supposing, for instance, the first and the last medium to 
be air, the emergent ray r’ will be, by (IV), 


r=r+2gN,_=P+ gy t+ ... +l, (22) 


where n, is the unit normal of «th plane face. In this more 
general case («>2) it will be convenient to draw all the normals 
N,, Ny, Ng and so on in one and the same sense, say, from the first 
towards the second, from the second towards the third medium; 
and so on. Again, as in the case of a single prism, the whole 
prismatic system can be represented by a skeleton of coinitial 
unit normals nj, Ng, Ns, etc., and unit edges e459, @y3, etc. These 
latter, no matter whether coinciding in space or not, can also be 
shifted to the same origin. (An air gap will also count as a com- 
ponent “prism,” with ~»=1.) Thus a system of N-1 prisms, 
i.e. of N interfaces, will be entirely represented by a skeleton of 
N vectors n, and N —1 edge vectors €;5, etc. If €, €s, etc., are the 
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prism angles, the relationship of the latter to the former unit 
vectors will be expressed by 
SIN €,.@j2=VnNg, sin €,. @93 = Vingng, etc. 

Any given incident ray r impinging upon this system will be 
converted into r’, as expressed by (22), where the g, are assumed 
to have already been found. (The system itself constitutes, of 
course, the best and more than complete framework of reference, 
thus giving no place for artificial ‘‘ auxiliaries.”) All properties of 
the emergent ray can at once be derived from the vector equation 
(22). We know already that the sole effect of each prism taken 
separately is to turn the incident ray round the edge as axis, 
leaving the angle 6 with the edge intact. The effect of the whole 
system will be the resultant of a number of such rotations round 
C12, €g3, etc., as axes, taken in their proper order, of course,—if, as 
in general will be the case, the edge vectors do not coincide, 7.e. if 
the prism edges are not all parallel to one another. 

The optical manufacturer attempts in ordinary circumstances 
(as in the case of spectroscopes) to make the prism edges parallel, 
4.€. €19=Co3=Cg4, etc. But, of course, he can succeed in doing so 
only approximately. Let us call such an ideal system a correct 
one, and any other, with edges more or less deviating from 
parallelism, a defective prism system.* Now, especially if a system 
of prisms is defective, formula (22) may be very helpful, because 
with n,, Ng, D3, etc., not coplanar, any of the ordinary Cartesian 
methods would soon lead to the greatest confusion and would 
require very careful drawings and possibly three-dimensional 
models to follow the rays in their tortuous course. The vector 
formula will therefore be particularly useful if it is required to 
estimate the effect of, say, small deviations from parallelism upon 
the emergent ray. 

If the system of prisms happens to be (or is considered as) a 
correct one, then all the normals n, are coplanar and all the unit 
edge vectors are fused into one, e, say. Multiply (22) by this 
vector ; then the result will be 

r’e=re, 
i.e. 6’=0, precisely as for a single prism, as could be expected. 
Whatever the orientation of the incident ray r, the prism-system 


* Although the reader may well imagine cases in which non-parallelism is a 
feature deliberately aimed at. 
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turns it rigidly round the common representative edge as axis, 
leaving its inclination @ intact. Notice also, that all the normals 
n, being coplanar, the vector sum in (22) is again a vector in their 
- plane. Such a vector can always be represented as the sum of 
two vectors along two fixed unit vectors in (or parallel to) that 
plane. Let these be a,b. Then (22) will assume the form 


r’=r+Aa+Bb, (23) 


where A, B are some scalar functions of the orientation of r, and 
of the characteristics of the prism system, 7.e. of the », and the 
relative orientation of the n, (the latter already including all the 
angles ¢,). These scalar functions A, B can easily be determined 
for the case of any concretely prescribed prism system. The 
resultant deviation (angle) ~=arccos(rr’) follows at once by 
multiplying (23) scalarly by r, which gives 
cos¥=1+A cos(r, a) +Bcos(r, b). (24) 
As a useful exercise the reader may treat fully the case of two 
glass prisms separated by an air-gap, either as an “‘ideal”’ or as 
a ‘ defective ” system. 


In the latter case it will be of particular interest to fully develop 
the formula for the effect of a small deviation from edge parallelism. 


7. System of Plane Mirrors. 


Let again, as in Section 6, each interface o, be plane, but let 
each of them be contemplated as a reflector. If the position of 
rays is disregarded, such a system, which may shortly be called 
a multiple reflector, will be entirely determined by the assemblage, 
and the order of succession, of the unit normals nj, ng, Ng, etc., 
drawn, let us say, away from the reflecting faces. 

As was already hinted at on p. 3, the resultant formula (IV) 
assumes for any such multiple reflector a particularly simple 
form. 

In fact, returning to the reflection formula (I), written down 
at the very beginning, and using now a single dash for every 
reflected ray, we have, for the first reflected ray, 


ry’ =1, —2n; (njr}), (25) 
where r, is the (original) incident ray. We already know that 
this is a particular case of the refraction formula (II), page 2, 
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viz. with p’ replaced by —y. But there is this important difference, 
that while (II) contained in its right-hand member, through g, 
the “ unknown ’’s’, formula (25) contains in its right-hand member 
r, only (besides the given n,) and no trace of r,’._ It is true that, 
as was promised on p. 2, the vector s’ can easily be eliminated 
from the right-hand member; but in doing so we obtain for the 
refracted ray a formula which in its structure differs essentially 
from the reflection formula (25). In fact, squaring (II), and 
remembering that s’?=p’2, s=y?, we have 


g°+2(sn)g=p?—-p?*, 
whence g, in terms of sn=4 cos7, and by substitution in (IJ), 

r’ =pr — pn (nr) +nV/o? +p? (ar)?, (26) 
where p=p/u’ and o*=1-p?. Herer only appears in the right- 
hand member ; but the operation to be performed upon r in order 
to obtain the refracted ray r’ is * manifestly a non-distributive one. 
That is to say, if (26) is written for the moment 

r’ =or; 
w standing for an operator, then w(A+B) is not equal wA +B. 
It is for this reason that we have said (Section 1) that it is not 
worth while to eliminate s’ from g. Nor would it give us any 
advantage to study the complicated and non-distributive operator 
w, as defined by (26),—which might be called the ‘‘ refractor.” 

On the other hand, a glance at the reflection formula (25) will 
suffice to show that the operator converting the incident into the 
reflected ray is distributive. This is the reason of the simplicity of 
multiple (plane) reflection, corresponding to a repeated application 
of such operators. 

In fact, n,(A+B) is n,A+n,B, which proves the distributivity 
of the operator 22,, say, involved in (25). That is to say, if we 
write that equation r;’={,r,, we have, for any A, B, 

0, (A +B) =0,A +2,B. 
This is a capital property of the ‘ reflector” 2, as contrasted with 
the non-distributive property of the “refractor” o. 

In order to see the structure of the ‘‘ reflector,” belonging to the 
mirror whose normal is n,, write (25), using the dot (instead of 
brackets) as separator, 

Tr,’ =F,—2N,.D4T,. 


* Quite apart from being more complicate than the reflecting operation. 
Cc 
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This is but a slight change in notation. Next, instead of writing 
twice the same operand r,, put 

ry =[1 —2ny.n]r,, 
to mean, by convention, precisely the same thing as before. Thus, 


the operator 2,, converting any incident ray into the corresponding 
reflected ray, is seen to be : 
0,=1-2n,.n,, (27) 
n, being the mirror’s unit normal. Such an operator is called a 
dyadic (n,.n, being called a dyad, since it contains fwo vectors 
separated by a dot). The reader need not, however, be deterred 
by new names or new concepts.* Since he knows how 2, was 
introduced, there can be no mystery about the meaning or the 
properties of this operator. In fact, he can read the definition (27) 
in plain English thus: to operate with 2, upon any given vector A 
means to cut off twice its projection upon n, along n, and to subtract 
the vector thus obtained from 1.A, that is, from A itself. This 
is a very simple operation indeed, more simple than most of those 
the reader daily performs, whether intellectually or physically. 
There will thus be no danger of a misunderstanding in using this 
operator, which is but a very particular example of the powerful 
dyadic. 
Thus our first reflected ray will be 

ry =0,7;. (28) 
Let now the reflected ray r,’ impinge upon the second mirror of the 
system. Then, 2,=1-2n,.n, being the “reflector” of this 
mirror, the second reflected ray will be, exactly as before, 


ry =Qar.=Qory’, 
and substituting r,’ from (28), 

Tr’ =2,2,7,, (28a) 
which means: operate with 2, upon r,, then with 2, upon the 
result of the first operation. This is in essence the same thing as 
to multiply a given number x first by 3 and then the result by 5. 
Both kinds are intelligible and reasonably simple operations. The 
only important difference between them is that while the order 
of the operations in the latter case is a matter of indifference, it 


*I mean, new to the reader. The name and the concepts, due to Gibbs, are 
old enough, 
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is not so in the former case. ,Q, is (in general) not the same 
operation as 2,95. So also is the actual reflection from mirror 1. 
followed by a reflection from mirror 2. not the same thing as a double 
reflection in the reversed order (unless the two mirrors are in a 
particular relation to one another). 

By (28a) it can now be said that the ‘reflector’? of a double 
mirror (nj, D4) is 

2,0). 

Similarly for three, four and any number of plane mirrors. 
The “ reflector’ of a multiple mirror, consisting of N plane mirrors, 
acting in the order I, 2, ... x, is 


Q =) OF; eee 2,00), 
or, more shortly, 


Q=T12,, where Q,=1-2n, .n,y, (29) 


the ‘‘ product” extending over all the component mirrors in their 
prescribed order. If any ray r impinges upon the multiple mirror, 
then the final ray, reflected in the said order, is 
r’ =Or. (30) 

The problem is thus reduced to building up the resultant operator 
Q from the component ones, 2,. Now, this is done almost as 
simply as the algebraic multiplication of a succession of poly- 
nomials. For to any succession of operators 2, belongs the 
associative property, 7.e. 

25(2,0,) = (2,0,)0), 

and similarly for any number of component reflectors.* The 
said simplicity of handling reflectors 2, and of building them 
up into the resultant reflector 2 for the multiple mirror is based 
upon these two properties, the distributivity and the associativity 
of these operators. In each case the process will consist in simply 
‘“‘ multiplying’ out the dyads n.n contained in the component 
reflectors. The only precaution needed is, to keep the order intact 
and to remember that juxtaposed vectors not separated by a dot 
are fused into ordinary scalar products. Thus, for example, 


N,. DNy. Ny=N,(NyNy) . Ny=AypNy « Dy, 
where a,,=n,n, =cos(n,, N.) isan ordinary scalarnumber. Similarly 


* Henceforth “reflector” will be used for the operator and ‘‘ mirror” for the 
reflecting surface to which it belongs. 
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for the product of any number of dyads. Such a “ product ” will 
always be again a dyad, including an ordinary scalar factor. 
In all the following developments we shall employ the short 
notation | 
DN, =COS (1, Nx) =Arey (31) 
so that, obviously, 
Qin =Anee 
Before passing on to examples, notice that equation (25) or (28), 
when squared, gives : 
1? =7,? +4(0yr,)? — 4(ayry)? =r. 
Thus, the incident r, being a unit vector, so is also the reflected ry’. 
The same property holds, of course, for any multiple mirror, 
i.e. if r’=Or, as in (30), then r'%=r?=1. Notice also that the 
iteration of any 2, is I or an idem factor ; in fact, 
0,0, =1 — An, . Dy +4(n, : n,)?=I — 4M, . Nye TDy . Dy 
(since n,2=1), whence 0,2, or 2,2=1, which is an obvious property. 
If y be the angle contained between the incident and the finally 
reflected ray, we have 
cos ¥=r(0r, (32) 
for any multiple mirror. The expression on the right hand is the 
scalar product of r into the vector Qr, or, which turns out to be 
precisely the same thing, of the vectorrQ intor. being a dyadic, 
composed of dyads a.b, can manifestly be used either as a post- 
factor ora prefactor. For any of the component dyads can be used 
so. Thus 
a.br means a(br) 
and ra.b means b(ar), 
both being perfectly definite results (vectors), in general differing 
from one another, of course. But r(Mr) is the same thing as (rQ)r, 
so that in (32) and all similar equations no brackets are needed. 
The reflector 2 belonging to any N-uple mirror will be obtained 
by multiplying out (with the said precautions) all the component 
reflectors {,. Thus, to obtain 2 for a quadruple mirror, say, we 
must necessarily have first the operators for a triple, and, before 
that, for a double mirror, Similarly for any N. It is therefore 
but natural to begin with the double mirror, and, before passing 
on to the triple and more complicated mirrors, to stop a while at 
each of the intermediate stages, discussing and interpreting the © 
results. It will be understood that by an ‘‘ N-uple” not neces- 
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sarily a system of N different simple mirrors is meant, but only that 
N successive reflections (in a prescribed order) are contemplated. 
An incident ray can be reflected from two mirrors, for instance, 
three or very many times, as is the case, for example, in the plan- 
parallel piece of a Lummer plate. 

Double Mirror. The cosine of the angle between the two com- 
ponent mirrors being a,.=nyn,, as in (31), we have, by (29), for 
the reflector of our double mirror, 

2=0,0, =[1 —2n, .n,] [1 —2n, . ny], 
1.é., with due regard to the prescriptions of the “ multiplication” 
rule, 
Q=1—-2n, .ny—2n,.n,+ 4ajn, .0j, (33) 
or, introducing the vector m =n, — 2a49N,, 
Q=1-2m .n,—2n,. Ny. (33a) 

This, of course, is valid for any operand, that is to say, for any 

incident ray r, giving r’ =Qr for the reflected ray. Thus 
r’ =F — 27,M — 274Ny, 
where 7,;=1N,, 7,=IM, are the given projections of the incident 
ray upon the mirror normals. Since the edge e, along which the 
two mirrors intersect, is perpendicular to n,, ny, and therefore also 
to m, we have again, as in the case of a prism, 
r’‘e=re=cos 0. 

In fact, the double mirror is but a sub-case of the prism (for p’ 
replaced by —y). The angle 6 which the ray makes with the edge 
remains intact after the double reflection. This is a general 
property of the double mirror independent of the value of a4p. 

From (33) we see that 2,2, differs from 2,, since the last 
dyad n,.n, is not symmetrical (ng, ny are two different directions). 
Thus, if a beam of parallel rays r is broad enough to impinge upon 
both mirrors, it is split by the double mirror into two beams, r’ as 
before and r” =2,0,r, so that 


r’ —r" =4a,,[Ny . Dy — Ny - NolT. (34) 
In particular, if the double mirror is orthogonal, that is to say, 
if dy,=n,n,=0, we have r’ =r", 


Z.e. the two reflected beams coincide. In this case the operator (33) 
becomes a symmetrical (or self-conjugate) dyadic, viz. 
2 =2,0, =0,0,=1 — 2n, . ny — 2M, . Ng. 
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Again, since n, 1 n,, we have 
n,.n,+n,.n,+e.e=I, 
for if i, j, kK, be any normal system of unit vectors (as is Ny, Ng, e) 
we have manifestly [i.i+j.j+k.k]V=V, for any vector V. 
Therefore, for an orthogonal double mirror, whose unit edge is e, 
Q=-—1+2e.e. (35) 
That is to say, the reflection from such a mirror is, independently 
of the order (12 or 21), wholly equivalent to the reflection from a 
simple mirror whose normal is e, followed by a simple reversal (— 1). 
This holds for any incident ray r. If the incident ray happens to 
be normal to the edge, 7z.e. if re=0o, we have 
r’=Or= -f 


? 

that is to say, the ray is sent back parallel to itself. This simple 
property is also familiar from ordinary geometrical constructions. 
But none of the cases of a skew r (oblique to e) will be found so 
easy to treat by the usual methods, based on a good drawing and 
spherical trigonometry, as by the proposed vectorial or dyadical 
method. 

Returning once more to the general formula (33), we have for the 
angle y between the incident and the finally reflected ray, by (32), 
on the other hand, since r=r,n, +7 n, + (re)e is a unit vector, 

142 +79 + 2ajo7 7. =1 — (re)? =sin? 6, 
so that the last equation can be written 
(I — cos ¥) =sin?@ =4a49747o, (36) 
giving the angle y=r, r’ for any double mirror and for any given 
incident ray r=7 ny +74No, 

From the general formula (36) we see without difficulty that there 
is no such double mirror which would send back (parallel to its 
own path) every incident ray ; in short, that a double mirror cannot 
be a ‘‘central”’ mirror, as are those used in modern signalling. 


In fact, cosy = —1 would imply 4a,97,7_ + (re)? =0, and this cannot 
be satisfied for all directions of r. 


Further discussions of the double mirror will be left to the reader 
as an exercise. In this case the adaptation of the vector formulae 


to numerical purposes is almost immediate, and calls, therefore, for 
no further explanations. 
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Triple Mirror. The resultant reflector in this case is 2 =02,0,0,, 
the product of 2;=1-2n,.n, into the operator of the preceding 
case already developed in (33). Thus, for any triple mirror, the 
order of reflections being 1, 2, 3, the resultant operator 2 =Q4og is 


Qieg =I = 2[n, . ny +N, . Ds +3 . Ns| 
+4[ayoMq . Ny +AggMg . Ng + Ag )Ng . Ny] — 8ay9%93Ng+ Ny. (37)* 


Here again, the second term being symmetrical, but not the third, 
nor the last, Q)59, etc., will differ from Q495. In other words, an 
incident beam of parallel rays r will give rise to six reflected beams 
represented by 
Diet, Oygor, Oogir, etc., 


which will, generally speaking, deviate from one another. These 
six reflected beams become parallel to one another, z.e. 2 becomes 
independent of the order of the three successive reflections when 


and only when 
412 = 493 = 43) =O,” 


abolishing the third and the fourth terms in (37), z.e. when the 
three component mirrors are mutually perpendicular. Moreover, 
Nj, Nz, Ng being in that case a triad of normal unit vectors, the 
dyadic n,.n,+n,.n,+n,.n3; becomes an idemfactor or I, and 
the operator is reduced to 


Gea —T 6... r =x. (37a) 


? 


That is to say, every incident ray, independently of the order of its 
three reflections, is sent back parallel to itself. The orthogonal 
triple mirror is a central mirror. Such triple mirrors, in fact, are 
used in practice, their orthogonality being now attainable to 
a high degree of approximation. 

Returning to the general triple mirror, write 7, 72, 73 for the 
direction cosines of the incident ray or its projections upon the 
three normals n,. Then the direction cosines of the reflected 
ray r’ will be linear homogeneous functions of 71, 7g, 7g, 012-, by (37), 


14! =[— 1 +4049? + dg” — 8449493493 |71 — [2442 + Ang g|%q — 244373, (38) 
and so on. 


*This may, at first sight, seem very complicate. Try, however, to treat the 
general triple mirror by the ordinary methods. Then (37) will be found very 


simple. 
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Exercises. Write down, by the general formula (37), the remaining 
two direction cosines ¥,’, 7,’, thus obtaining three scalar formulae 
with n,, N,, ns as reference system. 

Assuming that n,, n,, n, are not coplanar, 7.e. that the three 
mirrors constitute a proper pyramid (not a prism), any incident ray r 
can be put into the form 

r=P,N,+PN, +PsNs; 
where #,, etc., are some scalars. Determine them in terms of the 
orthogonal projections 7,=rn,, etc. (Notice that ~,=7,, etc., when 
and only when the mirrors are mutually perpendicular.) 

Similarly determine the coefficients of 

r’=p,'D, + paz + pss, 
and replace the formulae (38) by three equations giving the p,’ in 
terms of the p,, and in terms of the properties of the mirror, 7.€. as, 
etc. Discuss the general formulae thus obtained, and apply them to 
a number of particular cases. 

Treat the prismatic arrangement of the three component mirrors, 
when nj, N2, n3 are coplanar ; i.e. such that an, +bn,+cen,;=0, a, b, ¢ 
being three scalars, and so on. Consider first an incident ray | to 
the edges (all representable by one e), and then any oblique or skew 
ray. In the former case imagine the mirror pierced at a place, just 
to let through the incident ray; or else let the mirrors consist of 
non-silvered (plane parallel) glass plates. 


The angle ¥ between the incident and the finally reflected ray 
is given by cosW=rQr, 7.¢., according to (37), by 
B(I —cos f) = (ry? +79? +15?) — 2[ayar re + etc.] +4ay2e9%s1- (39) 
Each of the two bracketed expressions is symmetrical with respect 
to the suffixes I, 2, 3, but the last termisnot. Thus, if we replace 
the order of reflections 123, by 132, say, and call y’ the angle 


between the corresponding finally reflected ray and r, we have the 
elegant formula 


cos ¥ —cos f = 874A99(4497'3 — A451). (40) 
Similarly for any other permutation. Notice that 321 gives the 
same reflected ray as 123. 
If the three mirrors constitute a regular pyramid, i.e. if 
then formula (39) becomes 
(I — cos p) = (112 +142 +752) — 2 cos w(747g + 1a? +1973) 
+ 4775 cOStw, (39a) 
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More especially, if the incident ray r is equally inclined to the three 
reflecting planes of the regular aoe we have 


n=n="=3; cot - = 


and therefore, remembering that 1 — cos ¥ =2 sin? - 


sin¥ =cot 2/1 ~2 cos w +4 costa. (390) 
2 2 3 

In this case the angle y is, of course, independent of the order of 
reflections. The reflected beams (corresponding to a beam of 
parallel incident rays) are not parallel to one another, but are 
equally inclined to and symmetrically disposed around the direction ~ 
of the incident beam. These reflected beams will coincide in 
direction when and only when w = 90°, z.e. when the mirror becomes 
an orthogonal and, therefore, a central mirror. In fact, for 


o=5, formula (390) gives sin¥ =1, t.e. Y =180°, which means a 


reversal of the incident ray. This, as we already know, is the 
property of an orthogonal mirror with respect to any incident ray. 

Quadruple and more Complicated Mirrors. Multiplying (37) by 
) I—2n,. ny, where n, is the normal of the fourth mirror, the reader 
will obtain the operator Q2=Q,,.,, and similarly for quintuple 
reflections, and so on. As was already remarked, some of the 
normals n, may reappear, that. is, we may have, for instance, a 
quintuple reflection from but two mirrors; if such be the case, 
it is enough to put nz=n,;=n,, and ny=ny. 

As an example of this kind the reader may treat the rotation and 
progress of a skew ray incident upon an interior face of an equilateral 
prism built up of three plane mirrors; nj, nz, n; coplanar, ny=n,, 
Ns; =N,, and so on. ; 

Reversal of the Order of Reflections. Yor any multiple mirror, let 
r be the incident ray, r’ the finally reflected ray when the order of 
reflections is 123... N—1, N, and r” the finally reflected ray when 
the order of reflections is reversed, i.e. N, N—-1, .:. 321. Then, 
2 being 2,Qy_, ... 252,0,, as on page 25, we have | 

r=Qr and r”’=rQ, 
whence 
r’=rr, r’r=rr" =rfr, 
and therefore, for any multiple mirror, 
TY27 7; (41) 
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while rr” =19r, That is to say, the reflected rays 1’, r”, although 
not parallel to one another, are always equally inclined to the incident 
beam r. The equality ~42;=3q, mentioned on page 30 is but 
a special instance of this general property. 


8. Dyadic representing the most General Reflector. 


Consider the successive reflections of a ray r from any number 
of mirrors, whose normals are ny, ng, etc., and let 2, with the 
order of reflections subintended, be the corresponding resultant 
operator or the ‘reflector’? belonging to the whole, ordered, 
system. (This comprises also the case of any number of curved 
mirrors; the only difference being that when they are plane, all 
the normals are immediately given, and when curved, the successive 
normals are to be found by the application of the “ transfer for- 
mula,’’ developed in a previous Section.) Now 2 being the product 
of dyadics of the form 

I-—2n.n, 

is itself a dyadic. And since our space has but three dimensions, 
all the consequents * of any reflector, no matter how complicated 
the system to which it belongs, can be expressed-as linear homo- 
geneous functions of three non-coplanar vectors, which can always 
be made unzt vectors, let us say, f,g,h. The antecedents belonging 
to these consequents will then, for any given multiple mirror, be 
certain three vectors, let us say, a, b,c. Thus, the most general 
reflector will have the form of a trinomial dyadic, 


Q=a.f+b.g+e.h. (42) 


The plain meaning of this is that any incident ray r will be reflected, 
in the prescribed order, as 
r’ =a(fr) +b(gr) +¢(hr). (42a) 

The consequents f, g, h need not necessarily constitute a normal 
system; it may sometimes be convenient to choose an oblique 
system of vectors. The choice will depend upon the properties of 
the system of mirrors. 

The above result may also be stated by saying that the most 
general {2 is a linear vector operator, or r’ a linear homegenee\s 
vector function of r, as was already mentioned. 


* According to Gibbs’ terminology, the first vector, a, in a dyad a. b is called the 
antecedent, and the second vector, b, the consequent. 
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Splitting into Cartesians, for instance, and counting up the 
coefficients expressing the components of r’ in terms of those of r, 
the reader will see at once that the most general linear vector 
operator implies nine independent scalar or numerical data. Let 
us count the number of such data, implied in the most general 
reflector 2, upon the formula (42) itself. 

The consequents, f, g, h, which can always be taken as unit 
vectors, and, moreover, always as a normal system of unit vectors,* 
imply 2+1+0=3 scalar data (for the direction f being fixed, g 
requires but one datum, and h is codetermined, apart from the 
irrelevant sign or sense). Now it is always (i.e. for any trinomial 
dyadic) possible to give to the normal system f, g, h such an orienta- 
tion as to convert the antecedents also into some normal system of 
vectors,j but, of course, not unit vectors, after f, etc., have been 
made so. Thus the orientation of the system of antecedents 
implies three more data, and the size or absolute value of each 
antecedent one datum. Thus the total number of scalar data 
implied in the most general trinomial dyadic is nine. 

And since our reflector 2 is such that it leaves the size of any 
operand unaltered, 2.e. 

y2=—r7? 

for any r, the number of independent scalar coefficients contained 
in the most general reflector 2 is reduced to eight. These coefficients 
depend upon the properties of the system of mirrors and, of course, 
of the prescribed order of successive reflections. No matter, 
however, what their arrangement and how many of simple mirrors 
are built up into the system, all of its properties are crowded into 
eight coefficients,—corresponding to each prescribed order. On 
permutating the order, the operator, or its coefficients, will, in 
general, be changed; in some special cases they may remain 
intact. The reversal of order has, however, the perfectly general 
property expressed by (41). 

A dyadic, or linear vector operator, such that Qr=rQ, for any r, 
is called a self-conjugate or symmetrical one. The simple reflector 

* Although this may sometimes be not the most convenient choice. This, how- 
ever, does not influence the result of our counting. 

+A simple proof is given in J. W. Gibbs’ lectures on Vector Analysis, edited by 

“E. B. Wilson, New York, 1907, p. 304. Our dot, a separator, is not to be con- 


founded with Gibbs’ dot, which is the symbol of scalar multiplication. The 
notation adopted in the present volume is essentially that of Oliver Heaviside. 
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I—2n,n is, obviously, such an operator. But that corresponding 
to the general double mirror is no longer self-conjugate. Other 
details of the discussion of the multiple reflection formulae may 
be left to the reader. As to further examples and exercises, 
besides those proposed, the reader will doubtless invent them himself 
without difficulty. 


9. Differential Properties of Refracted or Reflected Pencils. 


Although the chief object of this book, viz. the tracing of a single 
ray, is in its essence exhausted by what has been said in the pre- 
ceding sections, it may be well to add one more section in order 
to give the reader at least some cursory hints about the treatment 
of differential properties of beams or pencils of rays sent through an 
optical system: properties, that is to say, which concern the 
passage from one ray of a pencil to its neighbours, or the passage 
from a given wave-length, and the corresponding index p, to an 
infinitesimally different one (astigmatism, chromatic aberration). 

Problems of this kind do not, as a matter of fact, require anything 
essentially different from what has already been given in the 
earlier part of the book. They are solved, in each case, by simply 
differentiating both sides of the vector formulae concernint a 
single ray, obtainable by the above method. Now, the differentia- 
tion of vector expressions does not differ from that of scalar ones, 
provided—of course—that we differentiate with respect to a 
scalar magnitude as the independent variable or parameter. As 
a precaution hardly anything more is needed than the preservation 
of order of the (vector) factors of a vector product and, in case 
of dyadics, the discrimination between antecedent and consequent. 

Thus, if r means the “‘ differential coefficient’ or the flux of r 
with respect to any scalar variable x, z.e. if r stands for | the 
flux of a sum a+b of two vectors is a +b, and that of the scalar 
product ab isab+ab. Similarly, for the vector product, 


is Vab = Vab + Vab, 
Ox 
and soon. Or, using “ differentials,”’ 
d(ja+b)=da+db, d(ab)=adb+bda, 
where da, db are infinitesimal vectors. Similarly 
dVab = Va db — Vb da, 
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the second term being negatived to compensate for the reversal of 
the order of its two factors. Further explanations, if at all needed, 
will be found in any treatise on Vector Analysis.* 

Such being the case, it will be enough briefly to consider the 
passage of a pencil of rays r through a single refracting surface 
(reflection being, formally, but a sub-case). Thus, returning to 
the fundamental refraction formula (II), Section 1, we have for 
any ray of an incident pencil, let us say for the ‘‘ central” ray of 
the pencil, 

s’=s+on, (11) 
where n is the normal of the refracting surface, of any shape (free 
of singularities), at the point struck by the central incident ray s. 
Let the refractive indices pw, »’ be fixed, and let us pass from the 
central ray s (to which corresponds the central ray s’ of the refracted 
pencil) to an infinitesimally neighbouring ray s+ds. If all the 
rays of the pencil are parallel, then, of course, ds =0, 8 representing 
the whole pencil. But in general we shall have ds#0. At any 
rate, since s*=y?, the relation sds=o, or dsts, will hold. The 
corresponding refracted ray will be s’+ds’, and ds’ 18’ as before. 
Now, passing from one ray to another means, in general, varying the 
normal n, and therefore also g. Thus, subtracting the original or 
central s’ from the new vector s’+ds’, we have for the required 
change of the refracted ray, 

ds’ =ds+gdn+n dg, (43) 
where, by the definition of g, p. 2, 

dg=psint, di—p'sin1’. dt’, 

and by the sine-law, »’ cos1’. di’ =. cos7. di, so that ultimately 
} dg = f(z) di, 
f being a known function of 7 alone.- Notice that, n being a unit 
vector, we have always ndn=o. In the case of reflection it: is 
enough to replace »’ by —p. 

In essence, formula (43) already contains the answer to all 
thinkable questions concerning the properties of the infinitesimal 
refracted pencil in terms of those of the incident one. It is enough 
to read it intelligently and to develop it appropriately in each 
particular case. 


* Cf, for instance, pp. 22-48 of my Vectortal Mechanics. Only a small part of 
what is there given will be required for the purposes of geometrical optics. 
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Let u, v be any curvilinear coordinates defining the position of 
a point at the refracting surface o. Then “‘ incident pencil given” 
means ds given as a function of u, v throughout an infinitesimal 
region round up, U, the point struck by the central ray. Now, the 
surface o being given (round that place at least), say, through the 
expression of its squared line-element, 


di? =E du? +2F dudv+G dv®, (44) 


dn is also known as a function of the two coordinates. The same 
is true of dg, which differs from di by a known factor f(7). 

Thus (43) finds ds’ as function of wu, v round up, vg; that is, it 
gives all the properties of the refracted pencil. For s’, n or 8’, 
Up, Vp give the direction and the position of the central pencil, and 

_ $/+ds’ with u+du, u)+du 
localize and determine the direction of any ray of the refracted 
pencil. Thus every question concerning the latter pencil is reduced 
to purely geometrical constructions. 


Employing u, v (for instance, ¢, 0 for a sphere) as our independent 
variables, we have, of course, 


Ot 
dg =f(t) (S dik +5dv), 
and in quite the same way 


on on 
dn =5,, 44 ae dv. 
If the incident pencil consists, for instance, of parallel rays, so 
that ds =o, then 
: ds’ =gdn +n dg, (45) 


whence nds’=dg, and Vnds’=gVndn, two convenient formulae 
replacing (45). The chief question in any such investigation is 
whether the neighbour ray passing through P(u,+du, v)+dv) 
meets at all, at finite or infinite distance, the central refracted ray 
passing through Po(u%, vo). In order to put this question in 
vectorial language, introduce the infinitesimal vector dl, drawn on 


the surface from Py to P, its size dl being determined by the quad-. 


ratic form (44) with E£, F, G as given functions of u, v. Then the 
said rays meet or not (are coplanar or not) according as the scalar 
magnitude 

Q =dl Vs'(s' + ds’) 
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vanishes or not. Remembering that Vs’s’=0, identically, this 
criterion-expression is reduced to 
Q=dl Vs'ds’, 
which, by (II) and (45), is 
Q=[dg Vsn+gVs dn+g?Vn dn|jdl. (46) 
Given the refracting surface and the incident pencil s, together 
with the central normal n, this expression will vanish or not accord- 
ing to the choice of the point P, 7.e. of the vectorial line element dl. 
In general we shall have Q =0 only for certain particular directions 
of dl upon the interface. Having established the conditions 
for Q=0, it will be a comparatively easy task to localize the point 
of intersection of the P-ray with the central ray of the refracted 
pencil, and thus to build up the caustics. If the incident beam 
is normal, then, whatever the shape of the refracting surface, 
the most convenient curvilinear coordinates u, v will be those 
measured along the lines of curvature of that surface drawn through 
Py as origin, in which case (44) is reduced to d/?=E du?+G dv. 
For oblique incidence, however, another network may in certain 
cases be found more convenient. 

As an exercise the reader may treat the case of normal incidence 
upon any refracting surface, such as an ellipsoid, a cylinder or 
a cone; and the case of oblique incidence for a spherical surface, 
thus deriving vectorially the well-known formulae for the sagittal 
and the meridional section of a pencil. 

Since these and similar problems have nothing to do with “ ray- 
tracing,’ the proper argument of this little volume, we cannot 
enter upon them any further. The last mentioned, and similar, 
exercises may, however, be warmly recommended to the reader 
as good and interesting opportunities of training in the use of 
vector language. 
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